Effects of coupling between octahedral tilting and polar modes on the phase diagram 

of PbZrx _3.Tij.O3 
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(Dated:) 

The results are presented of anelastic and dielectric spectroscopy measurements on large grain ce- 
ramic PbZri-zTLcOa (PZT) with compositions near the two morphotropic phase boundaries (MPBs) 
that the ferroelectric (FE) rhombohedral phase has with the Zr-rich antiferroelectric and Ti-rich FE 
tetragonal phases. These results are discussed together with similar data from previous series of 
samples, and reveal new features of the phase diagram of PZT, mainly connected with octahedral 
tilting and its coupling with the polar modes. Additional evidence is provided of what we interpret 
as the onset of the tilt instability, when is initially frustrated by lattice disorder, and the long range 
order is achieved at lower temperature. Its temperature Tit (x) prosecutes the long range tilt insta- 
bility line Tt [x) up to Tc, when Tt drops. It is proposed that the difficulty of seeing the expected 
i (111) modulations in diffraction experiments is due to the large correlation volume associated with 
that type of tilt fluctuations combined with strong lattice disorder. 

It is shown that the lines of the tilt instabilities tend to be attracted and merge with those of polar 
instabilities. Not only Tit bends toward Tc and then merges with it, but in our series of samples the 
temperature Tmpb of the dielectric and anelastic maxima at the rhombohedral/tetragonal MPB does 
not cross Tt, but deviates remaining parallel or possibly merging with Tt. These features, together 
with a similar one in (Na 1/ / 2 Bii/2) 1 _ a ,Ba I Ti03, are discussed in terms of cooperative coupling 
between tilt and FE instabilities, which may trigger a common phase transition. An analogy is 
found with recent simulations of the tilt and FE transitions in multiferroic BiFeOa [Kornev and 
Bellaiche, Phys. Rev. B 79, 100105 (2009)]. 

An abrupt change is found in the shape of the anelastic anomaly at Tt when x passes from 0.465 
to 0.48, possibly indicative of a rhombohedral/monoclinic boundary. 



I. INTRODUCTION 

The phase diagram of the most widely used ferroelec- 
tric perovskite PbZri-^Ti^Oa (PZT) still has unclear 
features (for the phase diagram, see Fig. [7]). It has been 
known since the fifties 1 - - — and the major recent discovery 
was the existence of a monoclinic (M) phased in a nar- 
row region at the morphotropic phase boundary (MPB) 
that separates the ferroelectric (FE) Zr-rich rhombohe- 
dral (R) region from the Ti-rich tetragonal (T) one. In 
the M phase the polarization may in principle continu- 
ously rotate between the directions in the T and R do- 
mains, so providing an additional justification for the well 
known and exploited maximum of the electromechanical 
coupling at the MPB. The existence of domains of M 
phase is actually still debated, the alternative being nan- 
otwinned R and/or T domains that over a mesoscopic 
scale appear as M££ Since experimental evidences for 
both types of structures exist, the possibility should be 
considered that genuine M domains and nanotwinning 
coexist at the MPB, being both manifestations of a free 
energy that becomes almost isotropic with respect to the 
polarization. 7 The part of the MPB line below room tem- 
perature has been investigated only after the discovery of 
the M phase, and is reported to go almost straight to K 
at x -0.52^ 

Recent studies are also revealing new features of how 



the TiC>6 and ZrOg octahedra tilt at low temperature. 
The instability of the octahedral network toward tilting 
is a common phenomenon in perovskites ABO3, usually 
well accounted for by the mismatch between the net- 
work of B-0 bonds with that of A-0 bonds which are 
softer and with larger thermal expansion! 14 ' 15 In these 
cases, lowering temperature or increasing the average B 
size sets the stiff B-0 network in compression, which is 
relieved by octahedral tiltingj 1 ^— In the case of PZT, 
Zr has a radius 19% smaller than Ti and one expects 
the tilt instability to occur below a Tt (x) line that en- 
closes the low-T and low-„ corner of the x — T phase 
diagram. Indeed, the Zr-rich antiferroelectric composi- 
tions are tilted (a~a~c° in Glazer's notation^ meaning 
rotations of the same angle about two pseudocubic axes 
in anti-phase along each of them and no rotation around 
the third axis), below a Taf (x) line that goes steeply 
toward K at x ~ 0.05. Also at higher Ti compositions 
tilting is observed {a~a~a~ compatible with the rhom- 
bohedral R3c structure) below a Tt line that presents a 
maximum at x ~ 0.16 and whose prosecution to low tem- 
perature was not followed beyond x — 0.4 until recently. 
The prediction from first principles calculations that oc- 
tahedral tilting occurs also in the M and T phased- has 
been confirmed by anelastic and dielectric^ structural^— 
Raman2_ and infra-red 2 ^ experiments. The presence of a 
low-temperature monoclinic Cc phased 2 , with tilt pattern 
a~a~c~ intermediate between tilted R and T has been 
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excluded by a recent neutron diffraction experiment on 
single crystals^ where below Tt coexistence was found 
of tilted R3c and untilted Cm phases. Yet, evidence for 
the Cc phase has been subsequently reported on PZT 
where 6% Pb was substituted with smaller Sr, in order 
to enhance tilting 

Here we report the results of anelastic and dielectric 
experiments at additional compositions with respect to 
our previous investigations, which reveal new features of 
the phase diagram of PZT, mainly related to octahedral 
tilting and its coupling with the polar degrees of freedom. 



II. EXPERIMENTAL 

Large grain (average sizes 15 — 30 /im) ceramic samples 
of PbZri-^Ti^Oa, with nominal compositions x = 0.05, 
0.062, 0.08, 0.12, 0.40, 0.487, 0.494 have been prepared 
with the mixed-oxide method in the same manner as pre- 
vious series of samples* 7 -^ The starting oxide powders 
were calcined at 800 °C for 4 hours, pressed into bars 
and sintered at 1250 °C for 2 h, packed with PbZr0 3 + 
5wt% excess Zr0 2 to prevent PbO loss during sintering. 
The powder X-ray diffraction did not reveal any trace 
of impurity phases and the densities were about 95% of 
the theoretical ones. The sintered blocks were cut into 
thin bars 4 cm long and 0.6 mm thick for the anelastic 
and dielectric experiments and discs with a diameter of 
13 mm and a thickness of 0.7 mm were also sintered only 
for the dielectric measurements. The faces were made 
conducting with Ag paste. 

The dielectric susceptibility x(cj,T) = %' ~~ was 
measured with a HP 4194 A impedance bridge with a 
four wire probe and an excitation of 0.5 V/mm, between 
0.1 and 100 kHz. The heating and cooling runs were 
made at 0.5 — 1.5 K/min between 100 and 800 K in a 
modified Linkam HFS600E-PB4 stage and up to 540 K 
in a Delta climatic chamber. 

The dynamic Young's modulus E(uj,T) = E' + 
iE" was measured between 100 and 770 K in vacuum by 
electrostatically exciting the nexural modes of the bars 
suspended on thin thermocouple wires^ The reciprocal 
of the Young's modulus, the compliance s = s' — is" — 
l/E, is the mechanical analogue of the dielectric suscep- 
tibility. During a same run the first three odd fiexural 
vibrations could be excited, whose frequencies are in the 
ratios 1 : 5.4 : 13.2. The angular frequency of the fun- 
damental resonance is^ w oc \f~E~ 1 \ and the temperature 
variation of the real part of the compliance is given by 
s (T) /so — Wq/w 2 (T), where luq is chosen so that so rep- 
resents the compliance in the paraelectric phase. The 
imaginary parts of the susceptibilities contribute to the 
losses, which are presented as Q^ 1 — s" / s' for the me- 
chanical case and tan S = x"/x' f° r the dielectric one. 



III. RESULTS 

For clarity, we will consider separately the anelastic 
and dielectric spectra with compositions in the range 
0.05 < x < 0.2, and those in the MPB region. We will 
present the new data together with those already pub- 
lished in Ref. GJ (x = 0.455, 0.465, 0.48, and 0.53) and 
Ref. (x = 0.1, 0.14, 0.17, 0.42, 0.45, 0.452). 



A. Octahedral tilting below Tt and Tit: 
0.062 < x < 0.2 

Figure [1] presents the dielectric and anelastic spectra 
measured during heating of PbZro.92Tio.osO3, a composi- 
tion where also the new transition at Tit is clearly visible 
in the elastic compliance s' . The comparison between the 
two types of susceptibilities puts in evidence their com- 
plementarity in studying combinations of polar and non 
polar modes. The dielectric susceptibility \' is of course 
dominated by the FE transition at Tq (note the loga- 
rithmic scale), it has a very attenuated step below the 
well known tilt transition at Tt, and practically noth- 
ing visible at Tit, due to both the broader shape of the 
anomaly and the proximity to the Curie- Weiss peak. The 
dielectric losses provide an indirect but more clear mark 
of the non polar transition at Tt, presumably through a 
change in the mobility and/or amplitude of charge and 
polar relaxations which arc affected by octahedral tilting. 

The effect of cooling through Tt on the dielectric sus- 
ceptibility is more convincingly shown to be a positive 
step in Fig. [21 and this fact will be discussed as a sign of 
cooperative coupling between tilt and polar modes. 

The elastic compliance s' , on the other hand, is only 
indirectly affected by the FE transition, since strain is 
not an order parameter of the transition and is linearly 
coupled to the square of the polarization. The Landau 
theory of phase transition ;] 27 ' 28 predicts a step in s' for 
this type of coupling, which is indeed observed at higher 
Ti compositions f^ 8 - but has a strong peaked component 
in Zr-rich PZT. We do not have an obvious explanation 
for this peaked response, which is frequency independent 
and intrinsic, but mechanisms involving dynamical fluc- 
tuations of the order parameter coupled with strain are 
possible^ The advantage of a reduced anelastic response 
to the FE instabilities is that the other transitions are not 
as masked as in the dielectric case, so that not only is the 
tilt transition at Tt clearly visible as a step in s' and peak 
in Q , but also the new transition can be detected at 
Tit even very close to Tq. As already discussed,— this 
transition has all the features of the transition at Tt but 
the associated anomaly is attenuated and broadened, so 
providing further support to an explanation in terms of a 
disordered precursor to the final long range tilt ordering 
below Tp. 

The broad peaks and steps in both the dielectric and 
anelastic losses below Tc have scarce reproducibility, 
which indicates their extrinsic origin, namely the motion 
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FIG. 1: Dielectric (left ordinates) and anelastic (right ordi- 
nates) spectra (higher panel real susceptibilities, lower panel 
losses) of PbZro.92Tio.08O3 measured during heating. 
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FIG. 2: (Color online) Dielectric susceptibility and loss of 
PbZri-^TirOs measured during cooling through the tilt tran- 
sition. 



of domain walls and charged defects, whose state depends 
on the thermal history. Instead, all the features indi- 
cated by arrows are completely independent of the mea- 
suring frequency, temperature rate and thermal history, 
and therefore are recognizable as intrinsic effects due to 
the FE and tilt transitions. Hysteresis between heating 
and cooling is observed due to the first order character 
of the transitions and to the presence of domain walls 
relaxations. Examples of the differences between the fea- 
tures that are intrinsic and stable and those that present 
dispersion in frequency or are less reproducible have been 
reported previouslyii^ and are omitted here. 

Figure [3] presents the anelastic spectra of PZT with 
0.062 < x < 0.17, including compositions already present 
in Ref. 7. All the curves are similar to the x = 0.08 case 
of Fig. [TJ with the three type of transitions at Tq , Tit 
and Tt clearly visible in separate temperature ranges. 
Both the s' and Q~ x curves have sharp peaks at the FE 
transitions, so that the Tc's are simply labeled with the 
compositions in %TL The other transition temperatures 
are indicated by vertical bars centered on the curves and 
arrows labeled with the respective compositions. The fea- 
tures of Q^ 1 in the Tit temperature range are not labeled 
because are due to the extrinsic contributions mentioned 
above. 

The temperatures of the tilt transition are identified 
with the upper edges of the steps in the real parts, which 
generally coincide with a spike or sharp kink in the losses. 



The rounded step and lack of reproducible anomaly in 
the losses increase the error on Tit, which however re- 
mains small enough to not change the features of the 
phase diagram discussed later. Due to the importance 
of the behavior of Tit (%) in the Discussion, a detail of 
this anomaly in the s' (T) curves, including 5% Ti, is 
shown in Fig 01 An anomaly corresponding to Tt might 
be present slightly above Tt for x = 0.17, but lacking a 
clear sign of it, it is assumed to coincide with Tt. At 
low x, the curve of x — 0.05 does not present any clear 
shoulder below To, and it is assumed Tt = Tq. 

The transition temperatures measured on both heating 
and cooling are reported in the phase diagram of Fig. 
where the Tt line departs from Tt at x ~ 0.17, has a 
kink centered at x = 0.11 and finally joins the Tq line at 
0.05 < x < 0.062. The new feature that will be the main 
focus of the present work is the kink and the merging with 
Tc at x > 0.05. It is also noticeable that the anomaly at 
Tt becomes more intense and sharper on approaching 



B. Compositions near the MPB 

As discussed in the previous investigations,^^ the 
MPB in PZT is signaled by a maximum in the dielec- 
tric and above all elastic susceptibilities. Again, it is 
stressed that such maxima are almost independent of fre- 
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FIG. 3: (Color online) Elastic compliance s' and energy loss 
coefficient Q" 1 measured at ~ 1.7 kHz on PZT at the compo- 
sitions 6.2, 8, 10, 12, 14, 17% Ti, as indicated by the numbers 
at the phase transitions. The curves of 10, 14 and 17% Ti are 
from Ref. Q. 

quency and temperature rate, and therefore are intrinsic 
effects due to the evolution of the order parameter at the 
MPB and its coupling with strain. Also the losses are 
rather high in the region of the MPB, but no feature is 
found that is directly ascribable to the phase transition; 
rather, their dependence on frequency and thermal his- 
tory show that they are due to the abundant twin walls 
and other domain boundaries, whose density and mobil- 
ity depend on many factors and is maximal around the 
MPB. Instead, the losses contain clear cusps or steps at 
the tilt transitions^^ so allowing Tt to be determined 
also in the proximity with the MPB, where the real part 
is dominated by the peak at Tmpb- We therefore discuss 
separately the real parts of \ and s, containing informa- 
tion on the polar transition at the MPB, and the losses, 
containing information on the tilt transitions. 



C. Maxima of the susceptibilities at the MPB 

Figure [5] is an overview of x' and s' curves measured 
during cooling at all the compositions x > 0.40 we tested 
so far. 

We call Tmpb the temperatures of the maxima in \' 
and s', marked with vertical bars on the curves. These 
temperatures do not coincide exactly with each other, 




FIG. 4: (Color online) Detail of the anomalies of the elastic 
compliance at Tit, measured at ~ 1.7 kHz during cooling on 
samples with 0.05 < x < 0.17. The numbers indicate the 
compositions in %Ti. The curves of 10, 14 and 17% Ti are 
from Ref. @. 

because x' an( i s ' are two different response functions 
of polarization and strain respectively, but, once plotted 
in the phase diagram, they present an excellent correla- 
tion with the MPB determined by diffraction, at least in 
the middle of the MPB line (see Fig. [7]). The dielec- 
tric maxima at Tmpb are much smaller and broader than 
the Curie- Weiss peak at Tc (note the logarithmic scale), 
whereas the anelastic maxima at Tmpb have comparable 
or even larger intensities than the step at Tq (part of 
the peaked component at To has a frequency dispersion 
denoting relaxation of walls «^). 

At x = 0.40 there is no peak attributable to the MPB, 
but only a minor step below Tt, which is indicated with 
triangles up to x = 0.452; beyond that composition, the 
step at Tt either disappears or is masked by the MPB 
peak. The other shallow anomaly centered at ~ 360 K in 
the curves up to x < 0.465 is the counterpart of the do- 
main wall relaxation appearing in the losses mentioned 
above and will be ignored. For x > 0.45 the peak at 
the MPB shifts to lower temperature and develops its 
maximum amplitude at 0.465, which has been argued to 
correspond to the point of the phase diagram where the 
anisotropy of the free energy is minimum^ The presence 
of a peak in s' at the MPB has also been argued to be 
evidence that the phase transition occurring at the MPB 
consists mainly in the rotation of the polarization, from 
the [001] direction of the T phase toward the [111] direc- 
tion of the R phase. In fact, in that case the transverse 
(perpendicular to the original [001] direction) component 
of P acts as order parameter and is almost linearly cou- 
pled to a shear strain, inducing a peaked response also in 
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FIG. 5: (Color online) Dielectric susceptibility and elastic 
compliance measured during cooling on PZT at the composi- 
tions indicated besides the curves in % Ti. The Tmpb's are 
indicated with vertical bars and Tr by triangles (only for 40 < 
x < 45.2). Present work: 40, 48.7 and 49.4% Ti; the other 
curves are from Refs. and ITU 



the elastic susceptibility, 7 ^ This would be an evidence 
that a monoclinic phase, and not only nanotwinned R 
and T phases, exists below the MPB. Yet, the smooth 
shape of the maximum is compatible with an inhomo- 
geneous M phase coexisting and possibly promoted by 
nanotwinning. 7 In fact, indications continue to accumu- 
late of intrinsic phase heterogeneity near the MPB com- 
positions also on single crystals^ 

Beyond x > 0.465, the peak at Tmpb gradually de- 
creases its amplitude and temperature, and, thanks to 
the great number of closely spaced compositions, is 
clearly recognizable as the signature of the MPB up to 
x = 0.515. The next composition, x = 0.53 (dashed 
curves), still has a shallow maximum at a temperature 
that prosecutes the Tmpb (x) line, but its nature appears 
different. In fact, the dielectric \' a t the same compo- 
sition lacks any sign of a maximum, and the overall s' 
curve does not any more prosecute the trend of the pre- 
ceding curves. For this reason, the temperature of this 
maximum at x = 0.53 is reported in the phase diagram as 
Tmpb but accompanied by a question mark. A Tmpb is 
extracted also from the curve at x = 0.42, even though a 
separate maximum is not present. It is however the only 



composition where s' has no sharp feature at Tq , and we 
assume that this is due to a rounded peak at Tmpb very 
close to To 



D. Tilt transition near the MPB 

The best signatures of the tilt transition below Tp are 
found in the anelastic losses. Figure [5] shows the Q^ 1 (T) 
curves at all the compositions x > 0.40 we tested so far 
(only 45.2%Ti is omitted in order to not overcrowd the 
figure) . 
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FIG. 6: (Color online) Elastic energy loss coefficient of PZT 
at the compositions indicated by the numbers (in %Ti), mea- 
sured during cooling at ~ 1.7 kHz. The anomalies at Tt 
are indicated by vertical bars and joined with a dashed line. 
Present work: 40, 48.7 and 49.4% Ti; the other curves are 
from Refs. Q and M- 

The Tp's up to x — 0.455 are the same as deduced 
from the step in the real part and indicated by triangles 
in Fig. [5] Up to x = 0.465 T T is identified with the 
temperature of a spike in Q (T), which gradually be- 
comes a cusp and starting from x = 0.48 becomes a large 
step. As in the previous figures, the Tp's are marked 
by vertical bars centered on the curves and joined by a 
dashed line, in order to better follow the evolution of 
the anomaly. The transition between the spike/cusp and 
the step anomaly is unexpectedly sudden, since it oc- 
curs within 0.465 < x < 0.48. Such a discontinuity ap- 
pears also in the dashed line joining the anomalies, and is 
marked by an arrow. We emphasize again that the losses 
generally have a limited reproducibility, because depend 
on the status of domain walls and defects; therefore, the 
regularity of the dashed curve joining the tilt anomalies 
of so many different samples is remarkable and testifies 
the good and uniform quality of the samples. 
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IV. DISCUSSION 

We refer to the phase diagram of PZT in Fig. [7J Be- 
low Tq and with decreasing Ti content, one finds the 
following phases : 13 ! 30 ' 31 ferroelectric (FE) tetragonal (T) 
PAmm with polarization P along [001], monoclinic (M) 
Cm with P rotated toward (111), rhombohedral (R) 
R3m with P || (111) and antiferroelectric (AFE) or- 
thorhombic (O) Pbam with staggered cations shifts along 
(110) and a"a _ c° tilt pattern. Below Tt octahedral tilt- 
ing occurs in all phases. 

In Fig. the solid lines join the transition tempera- 
tures deduced from our anelastic spectra measured dur- 
ing heating (filled triangles pointing upward), which are 
generally very close to the points deduced from the di- 
electric curves (empty triangles). The temperatures mea- 
sured during cooling are also shown as triangles pointing 
downward. The figur e contains all the data presented 
here and in Refs. I7lfl8l and, for completeness, also points 
obtained at compositions x < 0.05, that will be discussed 
in a future paper. 

The dashed lines are from the most widely published 
version of Jaffe, Cook and Jaffe^ with modifications of 
Noheda et aZ.— around the MPB. In a different version 3 ^ 
the Tmpb line below x — 0.45 does not prosecute straight 
up to join Tc almost perpendicularly, but rapidly de- 
creases its slope and joins Tc at x ~ 0.3. At present we 
have indication of a much smaller deviation of Tmpb from 
the datum at x = 0.42 (Fig. [5]), but already at x = 0.4 
there is no trace of a double transition. 



A. The octahedral tilt instability 

The tendency of the B06 octahedra in A 2+ B 4+ C>3 per- 
ovskites to tilt has been widely studied and can be ratio- 
nalized in terms of mismatch between a softer sublattice 
of longer A-0 bonds that compresses a stiffer sublattice 
of shorter B-0 bonds, until the incompressible octahedra 
tilt in order to accommodate a reduction of the cell and 
cuboctahedral volume Va without reducing their volume 
Vb- In the majority of cases the A-0 bonds are softer 
than the B-0 bonds, because they are longer and the 
cation A shares its valence with 12 nearest neighbors O 
atoms while B with only 6 of them 

This effect is quantitatively expressed in terms of the 
Goldschmidt tolerance factor 



t = 



A-0 



r\ + r Q 



V2lTO V2(r B 



ro) 



which is 1 for the cubic untilted case. In order to pre- 
dict a tendency to form a tilted phase, t is written in 
terms of the ideal ionic radii rx of the appropriate va- 
lence and coordination, which are tabulated^ 3 - If t < 1 
then the equilibrium B-0 bond lengths are too long to 
fit in a cubic frame of A-0 bonds, and tilting occurs 
below some threshold value; for example, many Sr/Ba 
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FIG. 7: (Color online) Phase diagram of PZT based on our 
anelastic and dielectric spectra. The solid lines join the 
anelastic data measured during heating; the dashed lines are 
those from Jaffe and Noheda. The question mark reminds 
that the shallow maximum of s' with x = 0.53 at that tem- 
perature probably does not signal the MPB crossing. 



based perovskites are tilted when t < 0.985 at room 
temperature^ Alternatively, the polyhedral volume ra- 
tio Va/Vb is defined^ which is 5 for the untilted case 
and becomes < 5 upon tilting. Of all the normal modes 
of a cubic perovskite, those that induce a decrease of 
Va/Vb are the combinations of rigid rotations of the 
octahedra about the cubic axes, while all other normal 
modes involve distortions of the polyhedra with little or 
no change in their volumes. 35 Rotations with all the oc- 
tahedra in phase along the axis are labeled M3, because 
the staggered shifts of the O atoms create a modulation 
with vector | (110), the M point of the Brillouin zone, 
while rotations with successive octahedra in antiphase 
along the axis are labeled R4, because the modulation 
is |(111), the R point of the Brillouin zone4£> These 
modes are also called antifcrrodistortive (AFD), because 
of the staggered modulation of the atomic displacements, 
and do not cause the formation of electric dipoles. Most 
of the low temperature structures of perovskites can be 
described in terms of combinations of these rotations to- 
gether with polar modes (shifts of the cations against the 
O anions). In particular, the low-temperature RZc rhom- 
bohedral phase of PZT, is obtained from the high temper- 
ature i?3m ferroelectric R phase, by applying anti-phase 
rotations of the same angle about all three pseudocubic 
axes (a~a~a~ in Glazer's notation^), while the O-AFE 
structure Pbam of PbZrOs is a combination of staggered 
and hence AFE displacements of the Pb 2+ and Zr 4+ ions 
along [110], with anti-phase tilts along [100] and [010] 
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(a-a-c ). 

In spite of the considerable amount of research on the 
phase diagram of PZT, little use has been done of these 
concepts in order to interpret it. An indication that the 
mechanism governing tilting in PZT is indeed the A-O/B- 
O mismatch is the observation that the substitution of 
6% Pb 2+ with the smaller Sr 2+ induces tilting in an origi- 
nally untilted T phase of PZT, while when codoping with 
Sr 2+ and Ba 2+ , the latter larger than Pb 2+ , the opposite 
effects of the two dopants on the average tolerance factor 
cancel with each other and no tilting is found. 37 When 
the tilt boundary Tt has been discussed in terms of t, 
the incongruence of the deep depression near the bor- 
der with the O-AFE phase has been pointed outp 3 - and 
explained in terms of frustration of AFE displacements 
of the Pb ions perpendicularly to the average FE direc- 
tion (111). Such displacements lack the long range order 
of the AFE-0 structure, and their frustration would be 
transmitted to the octahedral tilting through the Pb-0 
bonds, so lowering the Tt border in proximity with the 
AFE-0 phase.— Certainly the sharp depression in the 
border to the long range ordering of tilts, the Taf + Tt 
line, appears to be caused by some frustration, since it 
is a typical feature of the phase diagrams with compet- 
ing states , 38 ' 39 and both the FE/AFE polar modes and 
different tilt patterns may be in competition. 

Let us first consider the AFE/FE competition. The 
disordered displacements of Pb away from the average 
polarization have been proposed to occur over the whole 
R region of the phase diagram and particularly near 
the other MPB with the T phase, based on the large 
anisotropic displacement ellipsoids of Pb in Rietveld re- 
finements according to the R structure^ - or to the co- 
existence of M and R structures . 23 ' 30 While in these 
cases the Pb displacements off-axis with respect to the 
(111) direction may be imagined as having a FE cor- 
relation leading to a rotation of the polarization away 
from (111), recently a soft mode corresponding to AFE 
Pb displacements has been found near the MPB of PZT 
and in relaxor PMN and PZN-PT, suggesting that the 
AFE- like instability is a common feature of nanoscale do- 
main structures of rhombohedral or pseudorhombohedral 
lead-based perovskitesi 41 - Considering that such displace- 
ments occur over the whole R region or particularly at 
the MPBs with the O and the T phases, if they are so 
strongly coupled with the tilts, one would expect a de- 
pression of Tt also near the MPB with the T phase, which 
is not observed. There are some possible explanations for 
the different behavior of Tt at the two MPBs. One is that 
the polar displacements away from (111) might become 
static and with larger amplitude and AFE correlation 
only near the AFE border, so causing the frustration be- 
tween FE and AFE patterns, while at the MPB to the 
T-FE phase the AFE correlations remain at the stage of 
an anomalous phonon softening. 

On the other hand, it is possible that the main compe- 
tition occurs between the two different tilt patterns in the 
R and O phases, since the approaching of the tilt insta- 



bility to the FE one indicates that in the x — > region 
of the phase diagram polar and tilt modes have com- 
parable energetics. It is not evident, however, how the 
competition between a~a~c° and a~a~a~ tilt patterns 
would cause frustration, since they can transform into 
each other just by switching on and off the anti-phase 
tilt about the pseudocubic axis [001]. We believe that 
so strong a depression of the Taf + Tt lines is possible 
because partial tilting has already occurred at Tit- 



B. The intermediate tilt instability Tit (x) line 

Before discussing the nature of the instability at Tit, 
we emphasize the reasons why a positive step of the com- 
pliance during cooling like that at Tit should indicate 
a phase transition and not some kinetic effect related 
to domain walls or defects: i) cooling causes pinning or 
freezing of domain walls and therefore decreases the sus- 
ceptibility, while an increase is observed at Tt; ii) in the 
absence of tilting, the only conceivable walls just below 
Tc would be between the R-FE domains; if some anomaly 
in their behavior occurred around Tit, it would appear 
also, or mainly, in the dielectric susceptibility; Hi) the 
shape of the anomaly is independent of the temperature 
rate, history and frequency^ and therefore is an intrinsic 
lattice effect. 

The Tit line appears as the prosecution of the Tp one 
toward the lowest x, or equivalent ly lowest t, and highest 
temperature, and hence it has been identified as the on- 
set of precursor tilting/ Presumably, between Tit and Tt 
tilting would be disordered due to the enhanced disorder 
in the O-Pb-0 network near the AFE border, as sug- 
gested above.— The drastic depression of Tt would then 
be due to the fact that most of the mismatch between the 
Pb-0 and (Zr/Ti)-0 sublattices is relieved at Tit, and 
the final transition to a tilted phase with long range or- 
der requires the buildup of sufficient elastic energy in the 
disordered sublattice of tilted octahedra, that it is conve- 
nient to switch to the long range ordered phase. This is 
not alternative to the above discussion on the sharp min- 
imum of Taf + Tt in terms of competing phases. It is 
simply assumed that the frustration hinders tilting from 
reaching long range order but not from occurring on the 
local scale. 

This interpretation has not yet been corroborated by 
a structural study where the onset of tilting below Tit 
is actually observed and the tilt pattern is established. 
Actually, superlattice spots in electron diffraction have 
been interpreted in terms of in-phase A/3-type tilts within 
the otherwise Rim phase up to x ~ 0.15 i 42 ' 43 That re- 
gion roughly corresponds to the region delimited by the 
Tt border, though a border was not seen i 42 ' 43 These in- 
phase tilts could not be confirmed by x-ray or neutron 
diffraction and considerable debate ensued over the in- 
terpretation of the electron diffraction ^ (110) spots in 
terms of AFE (110) Pb displacements rather than Ms- 
type tilts, due to the much weaker strength of the latter 
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and the greater sensitivity of electron diffraction to the 
damaged surfacej 13 ' 44 




+ 

a a c Mo mode 



FIG. 8: Three tilt patterns of the BC>6 octahedra. The O ions 
at the vertices of the octahedra and the B ions at their centers 
are not shown; the A ions are white. The darker octahedra 
are in planes closer to the observer, a) a~a~a~ tilt pattern 
of the R3c structure seen from a direction close to the (111) 
rotation axis; b) out of phase and c) in phase rotation about 
the [001] direction. In b) and c) the octahedra of different 
planes share only O atoms that do not move under tilting.. 



The M3 tilts are in the first place in the search for 
a mechanism behind the phase transformation at Tit, 
because they would cause weak superlattice peaks at po- 
sitions coinciding with those from the disordered AFE 
Pb displacements, so explaining why their onset below 
Tit has not yet been noticed. In fact, the disorder in the 
Pb sublattice, including AFE-like shifts, survives even in 
the cubic phase, as demonstrated by the observation of 
electron diffraction | (110) spot s 42 ' 43 and the fact that 
EXAFS probes the same local environment of Pb 4 ^ and 
Zr/Ti^ both in the FE and in the cubic phases. It 
seems reasonable to assume that the AFE cation dis- 
placements develop their short range order together with 
the FE order below To, so producing superlattice peaks 
of the i (HO) type, that mask those from subsequent 
octahedral tilting below Tit, if this is of M3 type. The 
difficulty with the last assumption is that the final tilting 
below Tp is of R4 type, and the M3 in-phase tilting is an 
instability of different type rather than a precursor to it. 



C. § (110) reflections from hindered | (111) 
modulation of the tilts? 

In order to reconcile the observation of M-type mod- 
ulation below Tit with an expected short range i?-type 
modulation, we assume that the propagation of tilt fluc- 
tuations is hindered by the lattice disorder associated 
with Ti substituting Zr, with the consequent frustra- 
tion of the Pb displacements, and possibly other de- 
fects. Therefore, rather than to consider fluctuations of 
infinitely extended normal modes, it is more appropriate 
to consider the fluctuations of small clusters of octahe- 
dra, whose size is limited by the local disorder. This 
makes the different types of tilting represented in Fig. 
[5] inequivalent, since they have different correlation vol- 
umes. In a first approximation consider a network of 
rigid octahedra, that, in order to comply with the mis- 
match with the Pb-0 network, can only tilt without dis- 
tortions. This is the so-called rigid unit model, whose 
implications on the anisotropy of the phonon dispersions 
has been analyzed.— Here we focus on the effects that 
the anisotropy of the correlation length might have on 
the diffraction patterns. 

Consider first the rotation of a single octahedron about 
one of the cubic axes ([001] in Fig. [SJd) or c)); this will 
cause a rotation of all the other octahedra in the same 
plane perpendicular to the axis in an AFD fashion, cre- 
ating a \ (HO) modulation, with a correlation length l±, 
perpendicular to the rotation axis, which is an increas- 
ing function of the B-0 bond strength. In the pure rigid 
unit model l± is infinite and there is no correlation at all 
with the other planes of octahedra, because they share 
corners of the tilting octahedra only through immobile 
O atoms. In practice there is an interaction with the 
adjacent planes through the less energetic B-O-B bond 
bending and the extensions of the longer and weaker A- 
O bonds, resulting in a finite correlation length hi along 
the rotation axis, however shorter than l±. This means 
that the correlation volume surrounding an octahedron 
tilting about a (100) direction is a flat disc of diameter 
2l± and thickness 2iii with In <C and this will create 
superlattice reflections more intense and sharp at \ (110) 
than at | (111). On the other hand, the rotation of an 
octahedron about a (111) axis (Fig. Hk)), as the soft R4 
mode in the i?3c structure, will affect all the surround- 
ing octahedra, which share shifted O atoms both above 
and below the [111] plane containing the rotating octa- 
hedron. Therefore, for an octahedron rotating about a 
(111) axis the correlation volume is a sphere with diame- 
ter 2Zj_ , containing more octahedra than in the previous 
case. We therefore postulate that below Trr the magni- 
tudes of the tilts start becoming so large to cause their 
propagation through the correlation volume, as for any 
tilt instability, but this will occur first for those clusters 
where a deviation toward (100) tilt, with smaller corre- 
lation volume, is favored by the local cation disorder. In 
this first stage between Tit and Tp, only \ (110) super- 
lattice peaks would appear, easily masked by those from 
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AFE cation correlations. On further cooling below Tp, 
the increased instability of the R4 mode and long range 
elastic interactions trigger the long range tilt order of the 
R3c phase. In other words, in the absence of disorder or 
frustration Tp and Tit would coincide and represent the 
temperature below which the propagation of the a~a~a~ 
tilt instability starts. In the presence of disorder, local 
tilting about (100) axes is favored because of its smaller 
correlation volume; the resulting disordered tilting par- 
tially relieves the mismatch between Pb-0 and (Zr/Ti)-0 
bonds, and further cooling below Tp is necessary in order 
to reach the long range a~a~a~ ground state. This may 
explain the deep depression of Tx below x < 0.15. 



D. The kink in the tilt instability Tit (x) line and 
the effect of coupling of different modes on the 
phase diagram 

If octahedral tilts and polar modes were independent 
of each other, the Tc and Tx + Tit lines might approach 
and possibly cross each other in an independent manner. 
Instead, Tit merges with To at x — 0.06 with a noticeable 
kink around x ~ 0.1. The Tit line seems "attracted" by 
Tc , as if the tilt instability were favored by the ferroelec- 
tric one. On the other hand, also Tc may feel some effect 
from the proximity with Tit, since its slope slightly de- 
creases after joining with Tit below x ~ 0.05. It appears 
that both Tit and Tq increase with respect to the trend 
extrapolated from higher x, when they are far from each 
other. This is shown in Fig. where the dashed lines 
are the extrapolations (with no pretension of quantita- 
tive analysis) of Tq and Tp. The vertical arrows are the 
maximum deviations of the actual Tq and Tit lines from 
their extrapolations, AT C - 16 K and AT T ~ 48 K. A 
possible qualitative explanation of this observations is in 
terms of cooperation between a stronger FE instability 
and a weaker AFD tilt instability. The FE mode leaves 
the lattice unstable also below Tc, since its restiffcning 
is gradual, and also affects the modes coupled to it, in 
particular favoring the condensation of modes coopera- 
tively coupled to it at a temperature higher than in the 
normal stiff lattice away from Tq. The rotations of the 
octahedra are certainly coupled with the polar modes, as 
demonstrated by the polarizatio n 48 ' 49 and dielectrioZ^ 
anomalies at the tilt transitions, and the issue whether 
this coupling is cooperative or competitive is discussed in 
Sect. IIVHI Assuming that the first eventuality is true, 
one has a mechanism that enhances Tit in the proximity 
of Tc, but the coupling works also in the other sense: the 
proximity of the tilt instability favors the condensation of 
the FE mode, enhancing Tc. The energies involved in the 
FE instability are larger than those involved in tilting, as 
indicated by the fact that To > Tp, and accordingly the 
perturbation of the FE mode on the tilt mode is larger: 
ATt ~ 3 ATc . A phenomenological models that may be 
applied to explain these effects will be described in Sect. 
HTGl 
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FIG. 9: (Color online) Zr-rich region of the phase diagram 
based on our anelastic and dielectric spectra. The dashed 
lines are extrapolations from the high-x behavior of Tc (x) 
and Tr (x). The arrows ATc and ATr represent the deviation 
from the extrapolated behavior. 



E. Kinks in the Tmpb and Tr lines 

Other new features of the PZT phase diagram that 
derive from the present data are the approaches of the 
Tmpb (x) line with Tr and Tc , and a distinct kink in Tp 
when encounters Tmpb- This is better seen in the detail 
of the MPB region in Fig. [TUl where, besides the same 
data of Fig. [JJ other points of Tmpb and Tp are reported 
from the literature. The data are from diffraction 8 - (0) 
piezoelectric coefficient^ 1 - d\\ (□), 1/sn measured with 
piezoelectric resonance^ (o), Raman 9 - ( — ), dielectric (+) 
and infrared (x)2! spectroscopies. 

Let us first consider how Tp (x) enters the MPB region. 
The points from the literature, obtained from different 
techniques and samples, are rather sparse, but those from 
our anelastic and dielectric spectra have little dispersion, 
and show a clear change of slope of Tp (x) when it ap- 
proaches Tmpb at 0.487 < x < 0.494. This narrow com- 
position range is close to but not the same over which 
the anomaly in Q~ l changes between cusped and step- 
like (Fig. [6]). In fact that change, marked by arrows 
in Fig. [TUl occurs at x < 0.48, and therefore the two 
changes may depend on different mechanisms. We will 
discuss the transition in the shape of the Q^ 1 anomaly in 
Sect. IIVI[ and now we focus on the kink in the Tp line, 
which we think is closely connected with the proximity 
to the MPB. 

Also our Tvipb points draw a curve with little disper- 
sion, compared to the body of data in the literature, but 
in this case a difference emerges in the low tempera- 
ture region: even though with only three points below 
100 K, the data from diffraction 8 and Raman 9 define 
an almost straight MPB border that ends at T — at 
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FIG. 10: (Color online) Enlargement of the MPB region of 
the phase diagram of PZT. Triangles and solid lines from our 
measurements, as in Fig. [7] Dashed lines and diamonds from 
Noheda and other open symbols as indicated in the legend. 
The shape of the continuous Tmpb (x) line between x — 0.42 
and 0.45 is hypothetical. 



x = 0.520 ± 0.005. Instead, our data define a curved line 
that never crosses Xt- Our closely spaced points in the 
phase diagram and the regular evolution of the spectra 
from which they are obtained (Figs. [5] and |6|) suggest 
that the effect is real and characteristic of good quality 
ceramic samples. The last point with the question mark, 
obtained from the dashed curve in Fig. [5l probably does 
not correspond to Tmpb, but the difference remains at 
x = 0.515 between our curve at 120 K, and two points at 
50 — 60 K from diffraction and Raman scattering. These 
discrepancies may depend on differences in the samples 
rather than on the experimental technique. In fact, the 
existence of the intermediate monoclinic phase and its na- 
ture are not yet unanimously accepted, and it is also pro- 
posed that, besides nanoscale twinning, defect structures 
like planes of O vacancies may have a role in defining the 
microstructure of PZT and act as nuclei for intermediate 
phases* 5 ^. Hence, there is a range of microstructures that 
may well reflect in the position of the MPB, but, again, 
the consistency and regularity of the data encourage to 
consider the features presented here as intrinsic of the 
PZT phase diagram and not vagaries from uncontrolled 
defects. 

It results that also Xt and Tmpb almost coincide over 
an extended composition range, with Tmpb seemingly 
pushed up by Tp. For x > 0.49, Tmpb deduced from 
the maximum in s' and Tp deduced from the step in 
Q^ 1 run parallel and close to each other and it is dif- 
ficult to assess whether they still represent two distinct 



transitions or instead they are the manifestations of a 
same combined polar and tilt transition. 



Merging of tilt and polar instabilities also in 
NBT-BT 



Another example in which tilt instability lines merge 
with polar instability lines is (Nai^Bi^H-zBa^TiOa 
(NBT-BT). 54 This system has much stronger chemical 
disorder than PZT and a more complicated and less de- 
fined phase diagram, especially beyond the MPB com- 
position x (Ba) > 0.06, where the correlation lengths are 
so short to render the material almost a relaxor. In Fig. 
ITTI the phase diagram of NBT-BT is presented together 
with that of PZT. The broken lines represent the borders 
between regions with different tilts, while the solid lines 
are polar instabilities (Taf in PZT is both tilt and po- 
lar) . The regions were two types of instabilities merge are 
vertically hatched, while differently slanted hatches rep- 
resent different tilt patterns. In NBT the tolerance factor 
is small due to the smallness of the mean A ion size of 
Na + and Bi 3+ combined, and is increased by substitut- 
ing with Ba. Tilting occurs in two stages, first a°a°c + 
(T phase) below T\ and then a~a~a~ (R) below T2, and 
both Ti and T2 + Xvipb have negative dT/dx, so enclos- 
ing the low temperature/low tolerance factor region of 
the phase diagram, as discussed in Sect. IIV Al (the T\ 
line actually disappears into a highly disordered relaxor- 
like region). At variance with PZT, the polar instabilities 
occur at temperatures lower than those of tilting, and in 
two stages: first an almost AFE or ferrielectric region 
below a temperature T m signaled by a maximum in the 
dielectric susceptibility, and then FE below the so-called 
depolarization temperature T c \. The intermediate ~AFE 
structure is due to shifts of the Ti and Bi cation along 
[001] in opposite directions, so to make P almost null, 
and therefore is not the result of coupling with the AFD 
a°a°b + tilting.— The FE phase, however, appears below 
a Td < T<x at x = 0, which rises and merges with the 
decreasing T 2 (x) at x = 0.02 54 or 0.03^. Beyond 6% 
Ba, the Ti and T2 = Tmpb, the latter ill defined, sepa- 
rate again. As a result, the border Tj to the FE phase, 
instead of simply crossing T2, merges with it in the range 
0.02 < x < 0.06, following a wavy path. This would 
not be the only instance in the NBT-BT phase diagram 
where the tilt modes trigger a mixed tilt-polar transition, 
since it has been recognized that in pure NBT both the 
transitions at Ti and T2 are of such a type^ 

Similarly to the cases of Ti T /T c and T C /T MPB in PZT, 
the coincidence of Ti and T2 in the composition range 
0.02 < x (Ba) < 0.06 may be interpreted as a manifesta- 
tion of cooperative coupling between tilt and FE insta- 
bilities. Interestingly, like Pb 2+ also Bi 3+ has a lone pair 
electronic configuration, with the tendency to reduce its 
coordination number and form short bonds with covalent 
characte r 55 ' 58 which couple tilt and polar modest 
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FIG. 11: (Color online) Phase diagrams of PZT and NBT- 
BT,— where the broken lines are tilt instability borders, while 
the solid lines are polar instability borders. Vertical hatch- 
ing evidences the regions where the two types of lines merge. 
Differently slanted patterns represent different tilt patterns. 
P, F and A stand for paraelectric, ferroelectric and antiferro- 
electric; SRO and LRO stand for short/long range order. 

G. Merging of tilt and polar instabilities seen as 
trigger type transitions 

A possible mechanism for the merging of two transi- 
tions with order parameters (OPs) of different symme- 
tries had been proposed by rlolakowsky^ and carried on 
by Ishibashi;— This is also the case of octahedral tilting, 
whose OP is the rotation angle u>, and polar or antipolar 
modes with OP P (all one-dimensional for simplicity). 
Holakowski^ treated the case of a FE transition trig- 
gered by another OP, and we will adapt his arguments 
to a tilt transition triggered by the FE instability, refer- 
ring to Tit- In this case the minimal Landau expansion 
of the free energy is 

F = ?-p2+ b -P* + ^+ b ^+ C ^ + F c (1) 



where a = a(T — Tq) and a w = a u (T — TJ) represent 
the soft modes ideally vanishing at Tq and T u < To 
In the argument, the temperature dependence of a u is 
irrelevant, and Holakowski sets it constant; in our case it 
should be < Tit). The coupling part F c contains mixed 
terms P m oj n . Since u> and P have different symmetries, 
not all the mixed terms are invariant under the allowed 
symmetry operations in the cubic phase, under which 
F must be invariant; the lowest order term allowed by 
symmetry is^ 

F c = -\P 2 ^ 2 . (2) 

Such a term is always possible but generally overlooked 
when mixed terms of lower order prevail. If 7 > 
this term lowers the free energy when both P 7^ 
and describing a cooperative polarization-tilting 

coupling. After solving the equilibrium condition = 
dF/dP, the equilibrium value of u> is expressed in terms 
of equilibrium P, so that F is written only in terms of P, 
and F c renormalizes b as b' = b — j 2 /b^. The reason why 
the tilting free energy cannot be truncated to w 4 is that 
the renormalized b' u can become negative, in which case 
a positive 6th order term is needed to stabilize the free 
energy. We refer to the paper of Holakowsky for the de- 
tails and only report the result adapted to our case. The 
occurrence of the FE transition promotes tilting through 
the biquadratic coupling term, resulting in a tilt transi- 
tion at a temperature T t that can be also considerably 

higher than T u (which is zero in Ref. [60h . Increasing 
the magnitude of the coupling constant 7, the onset of 
the tilt instability is shifted to higher temperature and 
the following cases are encountered. When b' > the tilt 
transition is second order and occurs at 



hence a temperature higher than T w ,if T^, <^Tq. If 60 < 
b' < with bo = —4^/ then the transition is first 
order and occurs between T t and Tc; if b' < bo the free 
energy may have minima at both P / and uj 7^ 
(distinct from those of the pure FE phase with oj = 0) 
already at Tq and therefore a first order transition to a 
combined tilt/polar transition is possible at T > Tq. 

The above mechanism only requires that the tilt- 
polarization coupling is cooperative and enough strong, 
which is not forbidden by any constraint of symmetry or 
general principle, and we think that it can be at the basis 
of the anomalous rise and merging in temperature of the 
tilt and FE instabilities in Zr-rich PZT. The situation 
should actually be more complicated, since the triggered 
transition does not produce a phase with a clear sym- 
metry of the polar and tilt OPs, apparently because the 
energy shifts from chemical disorder compete with the 
energies involved in the regular Landau expansion of a 
homogeneous crystal. This disorder would be responsi- 
ble for preventing the complete tilt transition with long 
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range order down to Tt <C Tit, and should somehow be 
included into in the Landau expansion. These consid- 
erations should prevent from applying the above simple 
formulas for deducing the magnitude of 7 from the up- 
ward shift of Tit- Yet, an independent indication that 7 
is indeed large and positive comes from the positive step 
in the dielectric susceptibility below Tt, as shown in Sect 
IIVH3I 

The possibility should be explored that a similar mech- 
anism accounts for the proximity, instead of crossing, of 
the Tmpb and Tt lines at x > 0.50, as we observe in 
our large grain ceramic samples. In this case, the FE 
OP might be the transversal component P tl responsible 
for the rotation of P away from the tetragonal axis, and 
which takes a role in the peak of s' . 

The mechanism of the trigger-type transition has been 
applied so far to very few cases, like Bi4Ti30i2^ and 
recently proposed to explain the sequence of phase tran- 
sitions of the multiferroic BiFe03, 64 and is therefore con- 
sidered as very rare^ The possibility that a trigger-type 
mechanism is also responsible for the particular features 
of the phase diagrams of PZT and NBT-BT suggests that 
it may be not so rare. 

Finally, Tmpb seems to join also Tc smoothly, although 
the upper end of the Tmpb line in the phase diagrams 
above is based on only one datum and largely hypothet- 
ical (but in line with the much more marked effect in 
Ref. l32j). This case is different from the previous ones, 
since the OP active below Tmpb is not independent from 
that active at Tc, both being the polarization, and a 
triggered phase transition as above would be meaning- 
less. Nonetheless, the manner in which Tvipb meets Tc 
deserves further investigation. 



H. Competition or cooperation of tilt and polar 
modes? 

The interpretation above contrasts with the 
widespread opinions that the coupling between tilt 
and polar modes is negligible or competitive, and 
therefore it is opportune to discuss the nature of the 
interaction between polar and tilt modes. 



1. Negligible tilt-polarization coupling 

In the context of the general trends of the perovskite 
phase diagrams at different compositions, the coupling 
between tilting and deformation or polar modes has been 
considered negligible in rhombohedral perovskites. 65 
This assumption was based on the observation that in 
a large number of rhombohedral perovskites, including 
FE ones, the rotation angle uj of the octahedra and the 
polyhedral volume ratio mentioned in Sect. IIV Al are re- 
lated by 

Va/Vb = 6 cos 2 a; - 1 , (3) 



which is valid in the absence of octahedral distortions. 
In other words, the tilt angles are independent of the 
distortions^ which necessarily occur in the FE phases. 
This fact, however, shows that the mismatch between the 
bond lengths or polyhedral volumes is accommodated al- 
most exclusively by the tilting mechanism, and this is un- 
derstandable, since the only two normal modes producing 
a large change of Va/Vb are in-phase and anti-phase rota- 
tions of the octahedra, while all the other modes yield dis- 
tortions with little change of Va/Vb ^ Therefore, while 
the validity of Eq. is a manifestation of the fact that 
Va/Vb depends almost exclusively on the tilting angle oj, 
it does not imply a lack of coupling between tilting and 
the other distortion and displacement modes, for exam- 
ple influencing the type of tilt correlations. 



2. Competitive tilt-polarization coupling 

The most widely accepted view is of a competition be- 
tween FE and AFD instabilities, based on their opposite 
behavior under pressure i^ 6 - - — Indeed, in titanates and 
other perovskites pressure suppresses FE and promotes 
tilting, and this can be simplistically understood in terms 
of a steric mechanism: FE requires more space for cation 
off-centering while the octahedral rotation is promoted 
by a compression of the lattice that amplifies the mis- 
match between B-0 and A-0 sublattices. These observa- 
tions, however, show that the two instabilities behave in 
opposite manner under pressure, but not necessarily that 
they compete against each other. It should also be noted 
that there are several exceptions to the " rule" that pres- 
sure favors tilting, for example perovskites with trivalent 
B = Al;22 GdfZi and there are even studies that assume 
that the general rule is just the opposite: pressure induces 
the more symmetric cubic structure. 72 These different be- 
haviors under pressure depend on the relative compress- 
ibilities of the AO12 and BOg polyhedra, which can be 
rationalized in terms of the bond valence sum concept ^ 
It appears that in the zirconates and titanates the BOq 
octahedra are stiffer than the AO12 cuboctahedra, which 
is the usual case when B has larger valence than A. 

Another strong and widely accepted 69 indica- 
tion of competitive AFD-FE interaction comes from 
simulations 74 on SrTiOs showing that the FE mode 
of cation displacements along (111) and the AFD 
anti-phase tilt modes compete against each other. The 
simulations are done on SrTi03, which actually is not 
ferroelectric due to quantum fluctuations, but their 
result is clear: neglecting quantum fluctuations, the 
temperature of the tilt instability is 25% higher if the FE 
mode is frozen, while Tc is 20% higher if the AFD mode 
is frozen. This means that the two modes tend to cancel 
each other rather than cooperate, and the microscopic 
mechanism is identified in the mutual anharmonic 
interaction. The behavior of Tc and Tt in PZT or Ti 
and T2 in NBT-NT, however, seems the opposite. It is 
possible that the nature of interaction between AFD and 
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FE modes is different when the A ion is Pb instead of Sr, 
since the difference is not only of ionic size but there is 
also a more strongly covalent character of the bond when 
Pb goes off centerj 55 i 58 The other difference between 
the SrTi03 and PZT case is that in the first only the z 
component of the i?4 mode becomes unstable, leading 
to a T structure, while in PZT all three components 
yield the R structure. It is therefore unlikely that the 
simulation of SrTiC>3 can be plainly generalized to Pb 
compounds. Rather, it may be more appropriate to 
refer to similar simulations 64 on multiferroic BiFe03, 
where also Bi 3+ has the stereoactive lone pair like Pb 2+ . 
Such simulations indicate that the sequence of FE and 
AFD transitions is indeed dominated by the trigger-type 
mechanism, and that the coupling between the two 
types of modes is strong and both competitive and 
cooperative, due to the fact that the order parameters 
are multicomponenti^ 



3. Cooperative tilt-polarization coupling 

That the tilts are coupled with polar modes is demon- 
strated by the fact that below Xt there is a positive step 
in the real part of the dielectric susceptibility^ ' 18 ' 43 as 
clearly shown in Fig. [3J and in the polarization^ 8 - If 
the modes competed against each other, then the onset 
of tilting should depress rather than enhance the polar- 
ization and its derivative, and this becomes clear when 
considering the biquadratic coupling term, Eq. ©, the 
leading term allowed by symmetry. Below Xt the equilib- 
rium tilt angle uj starts growing, and F c renormalizes the 
term of the free energy oc P 2 as a' = a —juJ 2 and hence 
the dielectric stiffness = d 2 F/dP 2 ~ a'. Therefore, 
below Tr a positive step is observed in x if 7 > and 
a negative one if 7 < 0. From Fig. [2] it appears 7 > 0, 
namely the coupling between u and P is cooperative. In 
principle it would be simple to estimate the magnitude 
of 7 from that of the step in \' a t > the magnitude of 
ZJ from diffraction and a — a (T — Tq) from the Curie- 
Weiss peak. However, as noted in Sect. IIVF1 the simple 
free energy ([1]) does not contain the effect of disorder 
that depresses Tr, and the relevance of 7 deduced in this 
manner would be questionable. 

Additional indications of cooperative tilt-polarization 
coupling are the fact that the application of an elec- 
tric field in PZT modified with Sn and Nb enhances 
Xt of few degrees^ and its prediction in BiFe03 from 
a first-principle simulation* 6 ^ The coupling between FE 
and AFD modes has also been discussed in relation with 
the appearance of new Rama n 76 ' 77 and infra-red 21 - modes 
below Xt- 

We think that the fact that tilt and polar instability 
lines merge over extended composition ranges instead of 
crossing each other are due to a cooperative coupling be- 
tween polar/antipolar and tilt modes. It is interesting 
to note that a simulation of the PZT phase diagram in- 
cluding tilt degrees of freedom has already been donej 1 ^ 



and the Tmpb line, though finally crosses Xt, presents a 
marked bend on approaching it, exactly as appears from 
our anelastic and dielectric experiments. 



I. Transition in the shape of the Q 1 anomaly at 
Tt: a possible sign of R/M border 

As already noted in Sect. IIVE1 the kink in the Xt line 
and the transition in the shape of the Q^ 1 anomaly (Fig. 
EJ) appear at slightly different compositions, suggesting 
that the latter may have a different origin from the prox- 
imity to the MPB and hence polarization-tilt coupling. 

If this were the case, the most obvious explanation for 
the change of the Q _1 anomaly would be the postulated 
border separating R and M phases . 8 ' 9 ' 17 The existence of 
this border is one of the yet unsettled issues on the phase 
diagram of PZT, since there are various diffraction stud- 
ies, also recent and on single crystals) 23 ' 30 whose Rietveld 
refinements strongly suggest that the R and M phases 
coexist at least down to x = 0.4, so excluding a definite 
phase border. In addition, according to the view that the 
M phase is actually a nanotwinned R or T phase^i 6 - this 
border would not exist. Therefore, a R/M border would 
be highly significant: it would imply the existence of a 
long range M phase. A puzzling feature of this border 
would be its verticality. In fact, a truly vertical phase 
boundary in the x — X phase diagram would be under- 
standable at a specific composition that allows a phase 
to be formed with commensurate cation order. This is 
certainly not the case of the postulated R/M border in 
PZT, where no Zr/Ti ordering has ever been observed, 
and anyway x ~ 0.47 is too far from the closest relevant 
composition x = 5. Indeed, the R/M boundary found by 
first principles based simulations is not vertical: it starts 
at a triple point with Xc and Xmpb at x\ = 0.463 and 
ends at X = and X2 = 0.476.— No experimental evi- 
dence exists so far of the crossing of such a border with 
change of temperature, and the change of the shape of 
the Q _1 anomaly between 0.465 and 0.48 is not a conclu- 
sive evidence of its existence, since it might be associated 
with a change of the character of the transition through 
polarization-tilt coupling near the MPB. Further experi- 
ments at more closely spaced compositions are necessary 
to ascertain this point. 



V. CONCLUSIONS 

Anelastic and dielectric measurements are reported 
at compositions of the phase diagram of PhZri-^Ti^Os 
near the two morphotropic phase boundaries (MPB) of 
the rhombohcdral phase with the tetragonal and the or- 
thorhombic phases. Several new features are found in 
both regions, and discussed in terms of octahedral tilt- 
ing and cooperative coupling between the tilt and po- 
lar/antipolar modes. 
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We confirm the recent discovery^ of a new phase tran- 
sition at a temperature Tjt that prosecutes the border 
Xt of the tilt instability up to the Curie temperature 
Tc, in the region where Tt drops and meets the bor- 
der with the orthorhombic antiferroelectric phase. The 
new phase is assumed to represent the initial stage of 
octahedral tilting, without long range order due to the 
enhanced cation disorder near the AFE border. Lacking 
evidence from diffraction for this intermediate tilt region, 
the rationale for the assumption that tilting is involved 
is discussed in terms of mismatch between the networks 
of Pb-0 and (Zr/Ti)-0 bonds, as usual for tilted pcr- 
ovskites. In addition, it is proposed that, due to the 
anisotropy of the correlation length of different types of 
tilts, the initial stage of tilting in the presence of dis- 
order involves flat clusters of octahedra rotating about 
(100) axes, so producing i (110) type modulations, even 
when the final long range modulation is of \ (HI) type. 

The Tit tilt instability line merges with the ferroelec- 
tric Tc with an evident step and both temperatures ap- 
pear enhanced with respect to the extrapolations from 
the region where they are far from each other. Also the 
Tr line presents a clear kink when it meets the MPB and, 
contrary to previous experiments, Tmpb is found to de- 
viate and go parallel or even merge with Tr, instead of 



crossing it. These observations of deviations and merg- 
ing of tilt and polar instability borders are compared to 
a similar example in NBT-BT, and explained in terms of 
strong and cooperative interaction between the polar and 
the tilt modes, which causes a trigger type transition. 
Since the prevalent opinions are that tilt-polarization 
coupling is competitive or negligible, and the trigger-type 
transitions are extremely rare, the various indications of 
polar-tilt coupling in PZT are reviewed and discussed. 

Another feature that is considered is a rather abrupt 
transition in the shape of the anomaly in the elastic losses 
at Tt. The anomaly is a peak or cusp for x < 0.465 
and a step for x > 0.48. The possibility is discussed that 
between these two compositions there is an actual border 
between rhombohedral and monoclinic phases. 
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